D 


DOC  FILE  COP^  . ADA0644 1 8 


r 


» 


ELECTRICAL  CHARACTERIZATION  OF  GaAs 
ANODIC  OXIDES  USING  ISOTHERMAL 
DIELECTRIC  RELAXATION  CURRENT 


THESIS 

AFIT/GE/EE/78-18  Martin  J.  Biancalana 

Capt  USAF 


'•V'  k 'i 


iX  % 


O A 

\> 


145 


AFIT/GE/EE/78-18; 


& 


JLECTRICAL  CHARACTERIZATION  OF  GaAs 
ANODIC  CHIDES  j^SING  ISOTHERMAL 
DIELECTRIC  JELAXATION  £URRENT  * 


! 


'f'&H  ~ £ 'tk,  -e  £ / £_>  ( 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 


Requirements  for  the  Degree  of 


Approved  for  public  release;  distribution  unlimited 


Preface 


This  thesis  involved  the  design,  assembly,  and  use  of 
a measuring  system  for  dielectric  films  on  GaAs.  This  par- 
ticular system  was  decided  upon  after  unsuccessful  attempts 
to  use  capacitance-voltage  and  conduction  techniques  to 
characterize  GaAs  MOS  sample  devices.  The  bulk  of  the 
effort  went  into  producing  simple  and  reliable  apparatus, 
operating  under  control  of  a programmable  calculator.  The 
software  was  designed  to  be  easily  modified,  and  the  pro- 
cedure involves  no  critical  adjustments. 

Thanks  are  due  to  Dr.  Schuermeyer,  Dr.  Hartnagel, 

Dr.  Bayraktaroglu,  and  Mr.  Blasingame  for  their  very  help- 
ful discussions.  Dave  Mays  and  his  colleagues  at  the  Avion- 
ics Laboratory  provided  a great  deal  of  assistance  in  pre- 
paring the  samples.  Dr.  Borky,  Dr.  Fontana,  and  Prof.  Lube- 
feld  deserve  credit  for  reading  and  commenting  on  the  manu- 
script that  eventually  became  this  thesis.  Lastly,  a 
special  thanks  to  my  wife,  Jan,  whose  patience  and  encour- 
agement greatly  aided  in  the  completion  of  this  work. 

M.  J.  Biancalana 
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Abstract 


The  theory  and  practical  limitations  of  the  IDRC 
technique  of  Mar  and  Simmons  are  discussed.  An  improved 
version  of  the  original  experiment  that  incorporates 
digital  processing  of  the  data  is  described.  Results  were 
obtained  for  Si  and  used  to  establish  confidence  in  the  new 
procedure.  Anodic  oxides  of  GaAs  were  then  studied,  and 
the  results  indicate  that  there  are  large  densities 
(2  x lO^cm'^eV*1)  of  both  fixed  and  mobile  trap  states  at 
and  near  the  interface.  The  mobile  trap  states  behave  as 
though  negatively  charged.  Limitations  of  the  theory  as 
applied  to  GaAs  are  discussed,  and  a complete  software 
listing  for  the  procedure  is  included. 


ELECTRICAL  CHARACTERIZATION  OF  GaAs  ANODIC 


OXIDES  USING  ISOTHERMAL  DIELECTRIC 
RELAXATION  CURRENT 

I . Introduction 

Significance 

This  research  is  significant  in  that  it  adapts  a meas- 
urement technique  useful  in  Si  MIS  devices  to  the  problem 
of  measuring  the  electrical  properties  of  GaAs  passivation 
layers.  The  technique  of  isothermal  dielectric  relaxation 
current  (IDRC)  measurement  is  shown  to  provide  useful 
results  in  GaAs  measurements  despite  the  complex  nature  of 
the  GaAs  MIS  system.  This  work  was  done  as  part  of  the 
GaAs  surface  passivation  effort  being  conducted  by  the  Air 
Force  Avionics  Laboratory  with  the  goal  of  developing  devices 
for  high  speed  logic  and  microwave  systems  td  meet  a vari- 
ety of  projected  AF  operational  requirements. 

Surface  passivation  is  the  process  of  forming  an  insu- 
lating, protective  layer  on  a semiconductor  substrate.  An 
important  property  of  this  layer  is  its  ability  to  neutral- 
ize the  high  density  of  carrier  trap  states  normally  pres- 
ent on  a semiconductor  surface.  In  Si  technology  the  most 
common  passivating  layer  is  thermally  grown  or  deposited 
Si02. 

The  passivation  layer  is  the  key  to  fabricating  MOS  or 
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MIS  devices  with  any  semiconductor.  In  these  devices  the 
passivating  layer  forms  the  gate  insulation  and  accounts 
for  the  high  input  impedance  characteristic  of  the  device. 
Its  properties  also  strongly  influence  threshold  voltage, 
temperature  stability,  and  other  crucial  device  parameters. 

Were  MIS  technology  available  for  use  with  GaAs,  the 
high  electron  mobility  of  this  material  would  result  in 
high  speed  devices.  As  discrete  structures  this  would  mean 
more  efficient  microwave  devices  than  are  possible  in  Si  as 
well  as  achieving  significant  advantages  over  presently 
used  Schottky  barrier  gate  GaAs  devices.  In  integrated 
form  GaAs  MIS  circuits  are  expected  to  achieve  gigabit 
speeds.  Since  GaAs  cam  be  made  into  lasers  and  detectors 
already,  a capability  to  fabricate  MIS  circuits  with  this 
material  could  lead  to  complete  high  speed  communications 
systems  on  a single  substrate. 

Background 

Unfortunately,  surface  passivation  of  GaAs  is  a much 
more  difficult  problem  than  passivation  of  Si.  The  problem 
is  that  GaAs  is  unstable  at  the  temperature  needed  to  grow 
thermal  oxides,  and  that  As205  has  a very  high  vapor  pres- 
sure above  560°C.  When  As205  evaporates  from  native  oxide 
0»a2o3  is  left.  This  is  a porous  crystalline  material  that 
is  unsuitable  for  passivation. 

Layers  of  other  materials  such  as  SijN4  have  been 
sputtered  or  evaporated  on  GaAs  to  act  as  a passivating 
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layer.  These  are  again  too  porous,  and  leave  a very  high 
density  of  surface  trap  states.  Attempts  to  improve  the 
surface  state  density  problem  by  annealing  are  thwarted  by 
the  thermal  instabilities  of  GaAs  and  AS2O5,  and  the  poros- 
ity of  the  layers. 

Recently  another  technique,  anodic  oxidation,  has  been 
shown  to  be  very  promising.  Work  by  H.  Hartnagel  (Ref  l) 
and  others  has  shown  that  dense  layers  of  native  oxide  can 
be  formed  at  room  temperatures. 

This  method  involves  placing  the  GaAs  surface  to  be 
passivated  as  the  anode  in  an  electrolytic  cell  with  an 
inert  cathode.  A current  is  started  through  the  cell  by 
either  a current  or  voltage  source.  Whichever  parameter, 
current  or  voltage,  is  controlled,  the  other  is  monitored 
to  give  an  extremely  accurate  indication  of  the  progress  of 
the  anodization.  When  sufficient  oxide  has  been  grown,  the 
substrate  is  removed  from  the  electrolyte.  Since  the  common 
solvents  used  in  the  electrolysis  bath  are  ethylene  glycol 
and  water  or  ammonium  pentaborate  for  non-aqueous  work, 
this  is  an  inherently  low  temperature  process.  The  result- 
ing oxides  may  be  annealed  at  temperatures  as  high  as  350°C 
without  deterioration. 

In  addition  to  native  oxides,  more  complex  layered 
structures  have  been  formed.  Aluminum  sputtered  onto  the 
GaAs  substrate  before  anodization  forms  a layer  of  AI2O3 
during  anodization.  The  aluminum  oxide  forms  first  since 
A1  is  more  chemically  active  than  either  Ga  or  As.  Once 


3 


the  layer  has  been  formed,  it  is  possible  with  some  elec- 
trolytes to  form  native  oxide  either  over  the  AI2O3  layer 
or  under  it.  Growing  the  A1203  has  the  advantage  of  pro- 
ducing a better  dielectric  and  allowing  a higher  safe 
anneal  temperature. 

Anodic  oxidation  is  a simple  and  repeatable  process. 
During  the  process  the  reaction  rate  is  controlled  by  the 
driving  variable,  either  a voltage  or  current  source.  The 
progress  of  the  reaction  can  be  ascertained  by  monitoring 
the  other  variable,  either  current  or  voltage  respectively. 
This  allows  real  time  detection  of  abnormalities  such  as 
contamination  of  the  electrolyte.  It  also  greatly  aids  in 
achieving  repeatable  results. 

Problem 

' This  thesis  project  addresses  the  problem  of  determin- 

ing the  best  anodization  procedure  to  use  by  evaluating 
film  properties  to  demonstrate  empirically  which  anodic 
oxide  (or  any  other  insulator  for  that  matter)  produces  the 
best  dielectric  for  device  structures.  This  reduces,  in 
large  degree,  to  a problem  of  measuring  the  surface  state 
density.  In  Si  the  test  used  to  determine  this  is  a 
capacitance-voltage  measurement.  Surface  state  density  is 
found  from  these  measurements  using  a theoretical  model  of 
the  structure.  Conductance  measurements  rely  on  a similar 
model.  These  techniques  fail  when  the  assumptions  inherent 
in  the  model  are  no  longer  valid,  such  as  when  the  surface 
state  density  is  very  high  or  the  dielectric  becomes 
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lossy  (Ref  2:131).  Both  of  these  exceptions  occur  with 
GaAs. 

Surface  state  measurement  is  complicated  further  by 
the  lack  of  detailed  knowledge  of  the  surface  chemistry  of 
GaAs,  especially  under  anodic  conditions.  Also,  the  state 
of  the  art  of  growing  bulk  GaAs  is  still  not  comparable  to 
that  of  Si.  The  result  is  a high  level  of  surface  defects 
in  the  basic  crystal  lattice.  This  is  serious  in  itself, 
but  there  is  the  additional  possibility  of  adsorbtion  of  a 
doping  impurity,  such  as  carbon,  onto  the  surface  before  or 
during  anodization.  Complicating  any  speculation  about  the 
condition  of  the  surface  is  the  fact  that  carbon  adsorbed 
on  the  surface  could  be  at  a concentration  an  order  of  mag- 
nitude lower  than  the  level  which  is  detectable  by  Auger 
analysis  and  yet  produce  a surface  layer  doping  of  around 
10*^  states  per  cm-*  for  a depth  of  a few  atomic  layers. 
Carbon  impurities  produce  p-type  GaAs,  i.e.  very  effective 
trap  states,  in  otherwise  n-type  GaAs.  Many  such  conceiv- 
able hypotheses  about  the  GaAs  surface  have  not  been  tested. 
What  is  known  is  that  the  surface  state  density  is  extremely 
high  in  GaAs  MIS  structures  and  is  hindering  progress  toward 
MIS  devices.  It  would  be  of  great  value  to  have  a measure 
of  these  states  and  their  position  in  the  band  gap.  Such  a 
measure  is  the  goal  of  this  research. 

Approach 

The  approach  taken  in  this  work  was  to  adopt  a measure- 
ment technique,  based  on  rather  basic  concepts,  which  has 
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been  shown  to  work  on  Si;  assemble  it  in  improved  form; 
demonstrate  that  the  apparatus  works  as  expected  by  measur- 
ing a Si  sample;  and  finally  use  the  method  to  explore  GaAs. 

The  method  chosen  was  that  of  isothermal  dielectric 
relaxation  current  measurement.  This  technique  was  pro- 
posed by  Simmons  and  Wei  (Ref  3),  and  results  were  published 
on  MNOS  devices  by  Mar  and  Simmons  (Ref  2).  The  theory  is 
based  on  fundamental  characteristics  of  semiconductors. 
Therefore,  there  should  in  principle  be  little  trouble  in 
applying  it  to  GaAs.  The  method  also  should  work  well  when 
a high  density  of  states  is  present  (Ref  2). 

This  method  was  originally  implemented  by  its  devel- 
opers with  completely  analog  apparatus.  In  the  present 
apparatus  the  data  are  handled  digitally,  producing  a sim- 
ple operation.  The  digital  data  processing  also  allows  for 
straightforward  quantization  of  error  in  terms  of  sample 
standard  deviation.  The  result  is  a useful  measurement 
tool  that  is  simple  to  operate  and  understand. 
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II.  Theory 

Theoretical  Consider.-vtions  for  Silicon 

The  theory  proposed  by  Simmons  and  Wei  (Ref  3)  pro- 
vides the  basis  for  this  measurement.  It  is  an  attractive 
theory  in  several  respects.  The  physics  involved  is  fairly 
simple,  and  the  measurement  results  directly  in  a plot  of 
density  of  states  vs.  energy. 

The  theory  uses  the  definition  of  the  Fermi  level  and 
a description  of  how  the  Fermi  level  moves  in  energy  as  the 
bias  on  the  sample  is  changed  to  derive  an  expression  for 
the  discharge  current  as  a function  of  time  and  denro  ty  of 
surface  states.  In  an  unbiased  Si  MOS  structure,  the  energy 
distance  diagram  is  as  shown  in  Fig.  2-1.  The  Fermi  level 
is  flat,  and  due  to  surface  states  (typically  10*°  or  10** 
states  per  cm  ) , some  depletion  is  evident.  When  the  sample 
is  placed  into  accumulation  during  the  experiment,  the  dia- 
gram changes  to  look  like  Fig.  2-2.  Here  the  Fermi  level 
in  the  semiconductor  is  displaced  from  that  in  the  metal 
and  the  bands  are  bent  as  a result  of  the  excess  electrons 
at  the  interface.  The  shaded  area  represents  the  surface 
states  that  are  now  filled.  The  assumption  is  that  all 
states  below  the  Fermi  level  are  filled  during  accumulation 
and  the  rest  are  empty.  When  the  bias  is  changed  to  pro- 
duce surface  depiction,  the  diagram  tends  toward  that  shown 
in  Fig.  2-3.  This  change  cannot  happen  instantly  due  to 
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Fig.  2-2  Silicon  KOS  accumulated  by  external  bias. 
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the  time  constants  of  the  states.  Thus,  during  the  tran- 
sient the  states  in  the  top  half  of  the  band  (n-type  mater- 
ial is  considered  here)  discharge  in  sequence  moving  the 
Fermi  level  toward  the  configuration  in  Fig.  2-3.  This 
discharge  forms  a current  in  the  external  circuit  that  can 
be  measured. 

During  the  transient  the  Fermi  level  (Ept)  is  given  by 

EFt  = Ec-kT(ln*t  + *365)  ( 1 ) 

where 

Ec  = Energy  at  the  bottom  of  the  conduction  band 
k = Boltzmann's  constant 
T = temperature  in  °K 

J m attempt  to  escape  frequency  (=  v<rnNc) 

v = thermal  velocity  of  electrons 

0n  = electron  captive  cross  section  for  interface 
states 

Nc  = effective  density  of  states  in  the  C.B. 
t = time  in  seconds. 

This  non-steady-state  Fermi  level  is  approximated  by: 

E = E -kTlnJt  (2) 

P 

where  Ep  is  the  approximated  transient  Fermi  level. 

From  conservation  of  charge  considerations  the  follow- 
ing expression  for  the  current  density  due  to  the  discharge 
is  derived: 
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J = fEc  NS(E)  en  exp(-ent)dE 

liy 


where 


q = electronic  charge 

Cj  s insulator  capaci  tivnce 

= depletion  layer  capacitance 
Ev  = energy  at  the  top  of  the  valence  band 
en  = electron  emission  probability 

N,.(E)  = surface  state  density  as  function  of  energy 


By  noticing  that  en  exp(-ent)  can  be  approximated  by 
a delta  function,  the  integral  of  Eq  (3)  can  be  closely 
approximated  and  solved  in  closed  form.  With  the  assumption 
that  Cj^C^  an  expression  for  the  current  can  be  derived: 


I(t)  = Ns(Ep) 


where  A is  the  area  of  the  gate. 

This  equation  can  be  rewritten  as: 


Ns<Ep>  = 


qAkT 


to  show  that  the  density  of  states  at  the  energy  E is 

r 

directly  proportional  to  the  current- time  product. 

With  Eq  (2)  to  transform  the  time  axis,  an  l(t)*t  vs. 
t plot  can  be  rescaled  to  a plot  of  density  of  states  vs. 
energy. 
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Practical  Considerations 


There  are  several  aspects  of  the  theory  that  did  not 
explicitly  appear  in  the  original  paper  that  definitely 
affect  its  use  as  a basis  for  laboratory  measurements.  The 
first  of  these  is  the  fact  that  the  range  of  any  given  iso- 
thermal measurement  is  limited.  If  at  a given  temperature 
the  measurement  could  be  taken  from  10-"*  second  to  100  sec- 
onds the  range  in  eV  would  be  only  .28  eV  ( =10*^Hz ,T=200°K) . 
This  means  that  more  than  one  isothermal  measurement  must  be 
taken  to  explore  the  upper  half  of  the  band  gap  of  silicon. 
However,  there  are  more  practical  constraints  to  take  into 
account . 

If  there  are  any  states  closer  to  the  conduction  band 
than  the  range  of  the  measurement  (.24  eV  in  this  case), 
these  states  will  discharge  very  quickly,  and  their  dis- 
charge current  can  be  approximated  by  that  of  a capacitor: 

l(t)  = I0exp(-t/r)  (6) 

where  is  the  effective  time  constant  of  the  measuring 
system. 

When  this  is  plotted  as  l*t  vs.  t there  is  always  a 
maximum  at  t=T.  Therefore,  the  actual  measurement  must  not 
start  until  a time  much  greater  than  t=r.  Again  considering 
that  the  discharge  can  be  approximated  by  Eq  (6),  after  a 
long  time  the  current  will  be  too  small  to  measure  and  all 
that  will  be  in  the  system  is  noise.  That  noise  will  be 
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multiplied  by  an  ever  increasing  t (time),  so  the  noise 
will  look  like  an  increasing  or  wildly  erratic  density  of 
states.  These  considerations  lead  to  the  establishment  of 
a more  practical  measurement  window  of  10  seconds  to 
25  seconds. 

To  explore  the  states  over  the  whole  upper  half  of  the 
band  gap,  it  is  necessary  to  vary  the  temperature.  In  doing 
this  each  measurement  is  isothermal,  but  at  a different  tem- 
perature. This  makes  temperature  the  variable  that  moves 
the  measurement  window  to  the  place  in  the  band  gap  to  be 
explored.  At  liquid  helium  temperature  the  measurement  can 
be  as  close  as  .01  eV  from  the  conduction  band.  The  other 
end  of  the  measurement  is  always  the  midpoint  of  the  band 
gap  where  holes  produce  a significant  change  in  the  current 
(in  p-type  material  it  would  be  electrons  that  upset  the 
measurement) . 

As  a result  of  these  considerations,  the  apparatus  was 
designed  to  measure  the  average  isothermal  dielectric  relax- 
ation current  at  any  temperature,  calculate  the  sample  stan- 
dard deviation  of  the  averaged  data  sample,  and  form  the  I-t 
product.  This  is  then  plotted  as  a series  of  horizontal 
lines  whose  lengths  represent  the  energy  range  of  the  meas- 
urement, and  whose  vertical  positions  represent  the  average 
density  of  states  in  that  range.  A vertical  line  through 
the  midpoint  of  each  horizontal  line  extends  one  standard 
deviation  above  and  below  each  horizontal  line. 
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Considerations  for  GaAs 


There  are  known  differences  between  GaAs  MOS  structures 
and  their  Si  counterparts.  It  can  be  expected  that  some  of 
these  differences  will  affect  the  IDRC  measurement. 

The  first  significant  difference  is  that  it  is  not 
possible  to  bring  native  oxide  MOS  structures  into  accumula- 
tion (Ref  4),  At  a minimum,  this  has  the  effect  of  changing 
the  scaling  between  time  and  energy,  Eq  (2).  Since  Eq  (2) 
was  derived  from  the  occupation  function  which  is  part  of 
the  integral  that  defines  the  current  density,  Eq  (3),  it 
would  seem  that  the  current  density  function  is  also 
affected.  This  is  not  necessarily  so.  Since  the  integral 
in  Eq  (3)  is  evaluated  by  approximating  the  occupation  func- 
tion by  a delta  function,  Eqs  (4),  (5)  will  be  valid  as  long 
as  this  approximation  is  justified.  However,  while  the  den- 
sity of  trap  states  will  remain  proportional  to  I • t , the 
energy  scale  will  be  in  error.  While  the  error  will,  at 
the  least,  produce  a simple  translation  of  the  results  in 
energy,  it  could  also  introduce  expansion  or  contraction  of 
the  energy  scale. 

Another  difference  is  that  it  is  possible  that  the 
Fermi  level  is  pinned.  That  is,  the  structure  is  not  only 
prevented  from  reaching  accumulation,  but  is  also  prevented 
from  reaching  deep  depletion  (Ref  5:2),  Since  a wide  energy 
state  structure  has  been  observed  using  the  IDRC  technique, 
if  this  effect  is  present,  it  does  not  hamper  the  measure- 
ment. 
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The  effects  just  considered  are  caused  by  the  large 
density  of  trap  states  near  the  interface.  The  IDKC 
experiment  is  designed  to  investigate  these  states.  In 
the  present  form  of  the  apparatus,  the  IDRC  data  can  be 
easily  processed  to  yield  information  on  the  validity  of 
the  energy  scale  by  calculating  the  change  in  surface  po- 
tential represented  by  the  discharge  current  observed  and 
comparing  that  value  to  the  actual  change  in  bias.  If  the 
scale  is  valid,  the  bias  change  should  be  equal  to  or  larger 
than  the  calculated  change  in  surface  potential. 

The  change  in  surface  potential  is  calculated  as  a 
function  of  energy  by: 


where: 


VS(E) 


dE 


^ox 


cr^o 

d 


(7) 

(0) 


VS(E)  = change  in  surface  potential  as  a function  of 
energy 

Eq  = initial  energy 

CQX  = oxide  capacitance 

c = relative  dielectric  constant 

eQ  = permittivity  of  free  space 

d = thickness  of  the  oxide. 

In  calculating  Cox,  er  was  assumed  to  be  7. 
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III.  Experiment 


Equipment 

The  experiment  is  divided  into  three  parts  by  function 
or  operation.  It  is  therefore  convenient  to  describe  it  in 
terms  of  these  parts:  the  analog  section,  the  digital  sec- 
tion, and  the  temperature  controller.  These  sections  are 
shown  in  Fig.  3-1. 

The  analog  section  consists  of  those  instruments  that 
directly  interface  with  the  sample.  An  HP  6200B  power 
supply  is  used  to  produce  bias  voltages  of  plus  and  minus 
ten  volts  with  respect  to  ground.  A trigger  unit,  fabri- 
cated from  standard  parts,  is  used  to  adjust  and  switch  the 
bias  applied  to  the  sample  (Fig.  3-2).  A Fluke  8100A  digi- 
tal volt  meter  (DVM)  is  used  to  monitor  the  bias  voltage 
while  it  is  being  adjusted.  This  DVM  is  then  removed  from 
the  circuit.  The  trigger  unit  switches  bias  voltages  upon 
receipt  of  a trigger  signal  from  the  digital  section.  The 
leads  of  the  sample  header  are  soldered  into  the  circuit 
with  the  gate  lead  connected  to  the  output  of  the  trigger 
unit  and  the  substrate  lead  connected  to  the  picoammeter . 

The  picoammeter  is  a Keithly  417  Higi  ->eed  Picoammeter, 
and  is  connected  in  series  between  the  sample  and  ground. 

This  meter  has  an  output  for  a chart  recorder  which  is  used 
to  interface  with  the  digital  section. 

j 
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Fig.  3-1  Measurement  system  block  diagram. 
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The  digital  section  consists  of  the  equipment  connected 
directly  to  the  HP-1B  system  bus  (see  Fig.  3-1).  The  bus  is 
controlled  by  an  HP  9825A  calculator.  An  HP  9871A  printer 
and  an  HP  9872A  plotter  provide  the  output  of  the  digital 
section.  Data  are  entered  through  an  HP  3437A  DVM.  This 
voltmeter  acts  as  an  analog  to  digital  converter  and  as  the 
source  of  the  trigger  signal  for  the  analog  section.  This 
voltmeter  is  programmed  and  its  output  read  via  the  bus. 

The  temperature  controller  section  consists  of  a 
Varian  V-4540  temperature  control  and  a separate  Chromel- 
Alumel  thermocouple.  The  temperature  control  works  with  a 
stream  of  gaseous  nitrogen.  The  stream  is  first  cooled  in 
a bath  of  liquid  nitrogen  then  reheated  by  a feedback  con- 
trolled resistance  heater  to  the  desired  temperature.  The 
stream  leaving  the  heater  tube  is  constrained  to  flow  around 
the  sample  by  a shielded  styrofoam  insulting  jacket.  Since 
the  sample  is  several  millimeters  from  the  end  of  the  heater 
tube,  its  temperature  is  different  than  that  of  the  gas 
stream  in  the  tube.  To  accurately  measure  the  sample  tem- 
perature, a Chromel-Alumel  thermocouple  is  attached  with 
tape  to  the  sample  header.  The  output  of  the  thermocouple 
is  read  by  an  HP  3455A  DVM  capable  of  1 microvolt  resolu- 
tion. The  thermocouple  was  calibrated  in  liquid  nitrogen, 
and  its  corrected  output  converted  to  degrees  Celsius  by 
standard  tables  (Ref  6). 
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Sample  Preparation 

Samples  used  in  this  study  were  prepared  by  mounting 
them  in  TO-5  headers  with  conductive  epoxy  (Able  Bond  36-2). 
A gold  wire  was  then  connected  between  the  gate  and  the  lead 
of  the  header  using  the  same  adhesive  material.  This  pro- 
cedure is  necessary  due  to  the  fact  that  the  normal  bonding 
procedures,  ultrasonic  or  thermal  compression,  cause  severe 
damage  to  the  sample.  To  verify  that  this  procedure  intro- 
duced no  anomalies,  the  silicon  reference  sample  was  also 
prepared  in  this  manner. 


Procedure 

The  software  developed  for  the  calculator  (see  Appen- 
dix A)  handles  the  data  recording  and  processing.  This 
leaves  a very  simple  procedure  for  the  operator.  The  volt- 
age is  adjusted  using  the  controls  on  the  trigger  unit,  then 
the  temperature  is  selected  using  the  temperature  controller, 
and  the  device  is  biased  into  accumulation.  After  the  tem- 
perature as  measured  by  the  thermocouple  and  the  current  as 
measured  by  the  picoammeter  are  stable,  the  zero  of  the 
picoammeter  is  checked  and  the  temperature  recorded.  The 
routine  LONG  READ  is  initiated  causing  the  voltage  step  to 
be  applied  to  the  sample  and  the  data  to  be  recorded.  A 
second  reading  is  obtained  from  the  thermocouple  and  the 
temperature  average  is  computed.  The  data  is  then  scaled 
and  assembled  into  an  array  by  the  calculator.  The  array 
is  composed  of  500  ordered  pairs  (x,y)  where  x is  the  time 
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in  seconds  from  t=0,  and  y is  the  current  in  Amperes  at 
time  x.  Next,  the  routine  PACK  is  called.  It  calculates 
the  average  and  standard  deviation  of  the  current  readings 
starting  from  x=10  and  ending  at  x=2 5 seconds.  PACK  then 
adds  the  temperature  in  °K  as  a third  parameter  and  stores 
this  information  in  array  B. 

The  above  process  is  repeated  until  the  temperature 
has  been  stepped  across  the  desired  range,  A data  output 
is  then  obtained  as  a plot  (via  routine  PLOT  B)  or  a list- 
ing (via  routine  PRINT  B).  The  data  may  also  be  filed 
(using  FILE  B)  on  tape  for  future  use.  The  output  of  PLOT  B 
is  the  desired  profile  if  the  device  behaves  according  to 
the  theory  of  Simmons  and  Wei  (Ref  3). 

The  digital  system  makes  it  possible  to  save  even  the 
raw  data  arrays  (using  routine  FILE).  This  is  useful  in 
allowing  some  data  to  be  saved  and  processed  in  varying  ways 
without  repeating  the  experiment. 
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IV.  Results 


i 


Silicon  Results 

Reference  measurements  were  made  on  a Si  sample  to 

provide  a measure  of  confidence  that  the  system  functioned 

as  planned.  The  sample  was  n-type  Si  with  76  nm  of  ther- 

2 

mally  grown  oxide  on  a (ill)  surface  and  a .01824  cm  Au 
dot  as  a gate  electrode.  Two  measurements  were  made  on 
this  sample  using  different  bias  changes.  The  first 
(Fig.  4-1)  used  a change  of  +5  to  -5  volts.  The  second 
(Fig.  4-2)  used  a change  of  +8  to  -7  volts.  As  predicted 
by  the  theory  there  was  no  appreciable  difference  in  these 
results,  and  they  closely  matched  the  published  data  for 
density  of  states  vs.  energy  (Ref  7:452).  The  data  marked 
with  an  asterisk  (*)  are  the  result  of  a procedure  used  to 
detect  dependence  of  the  sample  on  bias  history.  This  pro- 
cedure will  be  explained  when  GaAs  is  discussed. 

Figure  4-3  shows  the  results  of  the  calculation  of 
change  in  surface  potential  with  energy.  This  calculation 
is  used  as  a check  on  the  validity  of  the  energy  scale. 
Since  the  maximum  change  in  surface  potential  is  less  than 
the  bias  voltage  change,  the  energy  scale  for  Si  is  in  fact 
as  assumed  in  the  theory. 

Gallium  Arsinide  Results 

The  GaAs  samples  studied  were  provided  by  B. 
Bayraktaroglu.  These  samples  are  native  oxide  MOS  devices 
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Data  for  Si  using  a +5  to  -5  volt  bias  change. 


Fig.  k-Z  Data  for  Si  using  a +8  to  -7  volt  change 


Fig.  4-3  Calculated  change  in  surface  potential  for  Si  data. 


produced  by  anodic  oxidation#  The  substrate  is  N-type  GaAs 

with  a doping  of  10  impurities  per  cm  . The  oxide  layer 

o 

is  100  nm  thick,  and  the  gate  is  a .00051  cm  A1  dot.  Two 
versions  of  this  device  were  studied.  Sample  B3A  was  not 
annealed  and  sample  B4A  was  annealed  at  350°  C for  15  min- 
utes. The  results  from  the  basic  experiment  are  shown  in 
Figs.  4-4,5.  The  data  point  marked  with  the  asterisk  is 
again  the  result  of  the  procedure  used  to  detect  dependence 
of  the  sample  on  bias  history.  The  procedure  was  to  cool 
the  sample  to  the  lowest  temperature  to  be  used,  then  warm 
it  in  steps.  At  each  step  a measurement  was  made.  A final 
measurement  was  made  at  a lower  temperature.  If  there  is 
little  or  no  dependence  on  bias  history  the  last  point  will 
fall  on  the  curve  obtained  during  the  heating  cycle.  This 
is  the  case  for  the  Si  data.  In  the  case  of  sample  B4A 
however,  the  last  point  falls  well  away  from  the  rest  indi- 
cating a dependence  of  the  characteristics  of  the  sample  on 
the  bias  history.  Figures  4-6,7  depict  the  results  of  two 
more  measurements,  first  with  a larger  bias  change  (+2.5  to 
-3  volts),  then  with  the  original  +2  to  -2  volt  change. 
These  results  support  the  idea  that  the  bias  and  bias  his- 
tory affect  the  characteristics  of  this  sample. 

The  plots  of  change  in  surface  potential  vs.  energy 
given  in  Fig.  4-8  for  sample  B4A  show  clearly  that  the 
energy  scale  is  in  error.  They  indicate  a maximum  change 
in  surface  potential  greater  than  the  applied  bias  change. 
This  is  not  actually  possible,  but  the  calculations  can 
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Unannealed  GaAs  data. 


Annealed  GaAs  data  using  +2.5  to  -3  volt  bias  change. 
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be  expected  to  give  this  indication  if  the  energy  range 
actually  explored  was  smaller  than  the  range  assumed  in  the 
calculation. 

Two  more  measurements  were  made  to  clarify  this  depend- 
ence on  bias  history.  Figure  4-9  shows  the  results  of  bias- 
ing the  sample  B4A  at  -2  volts  for  more  than  24  hours  at 
room  temperature , then  cooling  it  for  the  measurement  while 
maintaining  a bias  of  -3  volts.  The  bias  change  for  the 
experiment  was  +3  to  -3  volts.  Figure  4-10  shows  the 
results  of  3 hours  under  +3  volts  bias  at  room  temperature 
and  the  subsequent  cooling  of  the  sample  for  the  measurement 
while  it  was  under  +3  volts  bias.  The  region  of  special 
interest  in  Figs.  4-9,  10  lies  between  .35  and  .5  eV.  In 
Fig.  4-9  the  negative  bias  has  caused  a broad  maximum  in  the 
density  of  states  to  occur  in  this  region.  In  contrast, 

Fig.  4-10  displays  a minimum  in  the  region  of  interest  due 
to  the  action  of  the  positive  bias. 

The  most  plausible  explanation  of  this  behavior  is  that 
there  exists  a large  density  of  mobile  charges  in  the  oxide 
which  migrate  to  or  away  from  the  semiconductor-oxide  inter- 
face under  the  influence  of  the  bias  field.  Moreover,  since 
these  charges  appear  to  leave  the  interface,  i.e.  become 
less  effective  in  pinning  the  surface  potential,  under  a 
positive  bias,  it  appears  that  the  majority  are  negatively 
charged.  This  contrasts  sharply  with  the  familiar  situation 
in  Si  where  interface  charges  are  predominantly  positive. 
Given  the  lack  of  stoichiometry  in  the  oxide  on  GaAs,  it 
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may  be  that  Ga~  or  oth-»r  ionic  species  introduced  during 
anodization  are  responsible  for  this  mobile  charge. 

in  all  of  the  density  plots  there  is  an  indication  of 
a large  density  of  fixed  states  near  the  mid-gap  energy 
level,  and  some  of  the  results  show  fixed  states  near  the 
conduction  band. 
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V.  Conclusions  and  Recommendations 


Conclusions 

This  experiment  is  an  improvement  over  the  version 
used  by  Mar  and  Simmons.  While  the  basic  theory  is  the 
same,  the  implementation,  especially  the  digital  data  pro- 
cessing, allows  greater  flexibility,  and  provides  the  abil- 
ity to  perform  sophisticated  data  analysis  very  easily. 

The  experiment,  as  it  is  now  implemented,  provides  a 
powerful  tool  with  which  to  study  GaAs  MIS  structures.  The 
digital  data  processing  provides  the  flexibility  to  easily 
vary  the  way  the  data  is  collected  by  changing  such  param- 
eters as  the  sample  rate.  More  importantly,  the  data  can 
be  analyzed  quickly  and  accurately  to  yield  more  informa- 
tion. The  error  estimate  represented  by  the  standard  devi- 
ation is  one  example  of  such  an  analysis.  The  very  useful 
plots  of  change  in  surface  potential  vs.  energy  provide 
another  example  of  this  capability. 

The  results  of  measurements  using  this  technique  lead 
to  several  conclusions  concerning  the  anodic  oxide  MIS 
device  represented  by  Sample  B4A. 

1,  There  exist  high  densities  of  fixed  trap  states 
both  near  the  conduction  band  and  near  mid-gap. 

2.  There  exist  large  numbers  of  mobile  charges  in  the 
oxide;  when  located  at  the  oxide  interface,  they  introduce 
trap  states  that  lie,  in  energy,  between  the  two  peaks  of 
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fixed  trap  state  density 


3.  Due  to  the  large  density  of  trap  states,  both 
mobile  and  fixed,  the  Fermi  level  is  not  able  to  move  under 
applied  bias  as  assumed  in  the  basic  theory. 

4.  The  mobile  trap  states  respond  as  if  negatively 
charged,  they  are  mobile  at  room  temperature,  and  they  lie, 
in  energy,  between  the  peaks  of  fixed  trap  density.  These 
traps,  due  to  both  their  concentration  and  their  mobility 
make  this  oxide  unsuitable  for  device  application,  since 
they  would  make  MOS  thresholds  too  large  and  unstable  for 
useful  device  characteristics. 

The  results  of  the  change  in  surface  potential  vs. 
energy  plots  indicate  that  the  energy  scale  as  plotted 
according  to  the  theory  is  in  error.  The  scale  indicates  a 
range  of  energy  that  is  larger  than  the  range  actually 
explored.  This  means  that  the  transient  occupation  func- 
tion of  the  theory  will  have  to  be  modified  to  produce  a 
correct  energy  scale  when  used  with  samples  that  have  a 
large  density  of  interface  and  oxide  trap  states.  It  can 
also  be  noted  that  as  the  quality  of  anodic  or  other  oxides 
improves  to  the  point  where  device  structures  become  prac- 
tical, the  difficulty  introduced  by  large  mobile  charge  den- 
sities will  correspondingly  decrease,  making  this  measure- 
ment technique  increasingly  accurate  and  useful. 

Recommendations 

Based  on  the  success  of  this  experimental  measurement 
it  is  recommended  that* 
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1.  The  experiment  be  used  to  further  study  GaAs  MOS 
structures , 

2.  The  experiment  be  improved  with  the  use  of  a 
temperature  controller  capable  of  cooling  to  a lower  tem- 
perature , 

3.  The  theory  be  improved  to  yield  a more  accurate 
conversion  from  time  to  energy  in  the  presence  of  many  trap 
states,  specifically  by  incorporating  a transient  occupa- 
tion function  which  more  accurately  reflects  the  actual 
mechanisms  present  in  GaAs  anodic  oxides. 
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Appendix  A 


Calculator  Software 

Introduction 

The  data  collection  and  processing  for  this  experiment 
was  done  with  a Hewlett-Packard  9825A  calculator.  The 
String  Variable-Advanced  Programming  (98210A),  Matrix 
(98211A),  and  Plotter-General  I/O-Extended  I/O  (98216A) 
read  only  memories  were  installed.  The  calculator  has  an 
integral  magnetic  tape  cartridge  to  store  program  and  data. 
Upon  inspecting  the  instruction  set  for  the  tape  cartridge, 
it  became  apparent  that  if  large  amounts  of  data  or  complex 
programs  were  to  be  used,  the  operation  of  the  tape  car- 
tridge system  would  have  to  be  as  automatic  as  possible.  A 
system  program  called  TAPE  MONITOR  was  written  to  do  this. 
This  system  program  and  the  program  for  the  experiment,  were 
written  in  modules  or  states.  This  allowed  the  easy  addi- 
tion or  deletion  of  program  modules  as  new  ideas  were  form- 
ulated and  old  ones  discarded. 

Variables  and  Conventions 

To  pass  data  between  the  modules  some  conventions  were 
established.  The  following  "r"  variables  were  given  special 
meaning} 

rl  = current  number  of  tape  files 
r2  = current  track 
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r3  = current  file 

r4  = data  array  length/byte  length  for  data  or  program 
r5  = status  variable 
r6  = scale  factor 
r7  = scale  factor 
r8  = temporary  storage 
r9  s automatic  mode  indicator 
rll  = the  next  usable  program  line  after  the 
experiment  program 

rl2  = the  next  usable  line  of  program  after 
TAPE  MONITOR 

(rll  and  rl2  are  used  to  link  program  modules  from 
tape  memory  to  calculator  main  memory) 
rl5,  rl6,  rl7,  rl8,  and  rl9  are  used  in  the  CONVERT 
subroutines . 

A string  variable  matrix  and  a numeric  matrix  are  used 
to  store  the  file  tables  needed  in  the  system  routines. 
F[50,4]  is  allocated  as  follows: 

F[nfl]  track  number 
F[n,2]  file  number 

F[n,3]  status  where  2 is  used  and  0 is  free 
(F[l,3]  is  used  to  store  rl ) 

F[n,4]  absolute  size  of  the  file 
F$[50,lo]  contains  the  names  of  the  files.  Each  name  can 

i'* 

be  up  to  lO  characters  long. 

String  variables  A$  and  B$  are  used  to  pass  names  and 
labels  respectively. 
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A[500,2]  is  an  array  of  ordered  data  pairs  500  pairs 
long.  B[40,3]  is  an  array  or  ordered  data  triples  40 
triples  long.  The  arrays  A and  B are  used  to  transfer  data 
from  one  module  to  another. 

Each  module  is  described  separately  in  as  much  detail 
as  its  complexity  warrants.  Verbal  descriptions  and  flow 
diagrams  are  used  as  needed  to  explain  the  listings. 

System  Program 

The  system  program,  TAPE  MONITOR,  is  really  more  of  a 
file  manager,  and  is  composed  of  modules  that  aid  in  using 
the  tape  cartridge  memory.  The  first  module  (Figs.  A-1,2) 
is  an  initialization  and  selection  or  directory  state.  It 
dimensions  the  arrays  needed  for  the  operation  of  the  sys- 
tem modules  in  a manual  mode.  Two  of  the  modules  are  vec- 
tors to  the  experiment  program,  DATA  MANIPULATOR.  These 
transfer  control  to  the  experiment  program  which  retains 
control  unless  it  is  directed  by  the  operator  to  return 
control  to  the  system.  The  other  modules  all  return  control 
to  the  directory  when  they  are  in  the  manual  mode.  This 
distinction  is  made  since  MARK  FILE  and  LOCATE  FILE  are 
capable  of  being  used  as  subroutines,  and  in  this  automatic 
mode  they  return  control  to  the  calling  routine. 

PRINT  TABLE  and  FREE  FILE.  The  simplest  of  the  system 
modules  are  PRINT  TABLE  and  FREE  FILE  (Fig.  A-3).  Each  of 
these  does  just  what  its  name  implies.  PRINT  TABLE  prints 
a listing  of  the  current  status  of  the  tape  file.  It  finds 
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J*'ig.  a-1  Tape  Monitor  flow  diagram. 
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print  table  routine 
: wr  t 7 02,  " ";v;rt  702,"  " 

2,5x , "NAME" , 5x , "TRK" , 5x , "FILE " , 5x , "Current  Size 


A-3  Listing  for  FRINT  TABLE  and  FREE  FILE  modules 


part  of  its  information  in  F and  F$,  and  finds  the  rest  by 
using  an  "identify  file"  command  in  the  calculator  language. 
FREE  FILE  asks  for  the  name  of  the  file  to  be  freed,  then 
checks  for  its  existence.  If  it  exists,  the  file  space  is 
marked  as  free,  and  its  name  is  erased.  If  no  such  file 
exists,  an  error  message  is  displayed,  and  control  is 
returned  to  the  directory. 

MARK  FILE  and  LOCATE  FILE.  The  next  two  states,  MARK 
FILE  and  LOCATE  FILE,  are  more  complex.  Each  has  two  main 
flow  paths.  One  path  is  used  for  the  manual  mode.  It 
interacts  with  the  operator.  The  other  path  is  an  automatic 
mode  that  operates  as  an  invisible  subroutine,  that  is, 
there  is  no  interaction  with  the  operator  to  identify  that 
it  is  running.  The  flow  diagrams  are  the  easiest  way  to 
trace  the  flow  through  these  routines.  In  the  automatic 
mode  A$  is  assumed  to  have  the  file  name,  and  r4  is  assumed 
to  have  the  proper  length.  In  the  automatic  mode  r9  must 
be  set  to  1. 

MARK  FILE  looks  for  a free  file  long  enough  to  accommo- 
date the  new  file.  If  the  space  is  not  found,  a new  file 
is  marked;  the  file  table  is  updated  and  recorded;  and  the 
track  and  file  numbers  are  passed  using  r 2 and  r3.  If  the 
space  is  found,  the  files  are  updated  and  recorded  and  the 
track  and  file  numbers  are  passed.  (See  Figures  A-4,5.) 

LOCATE  FILE  finds  the  track  and  file  numbers  of  any 
active  file.  If  it  is  asked  for  a non-existent  file,  an 
error  is  indicated.  When  the  file  is  found,  the  track  and 
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Fig.  A-6  Flow  diagram  for  LOCATE  FILE 


FILE 


file  numbers  are  passed  using  r2  and  r3  (See  Figures  A-6,7.) 


Experiment  Program 

The  experiment  program,  DATA  MANIPULATOR,  is  also 
written  in  modules  or  states.  Again  they  are  loosely  con- 
nected by  a directory  or  selection  state.  One  interesting 
feature  is  that  some  of  the  modules  allocate  memory  in  some 
way,  and  therefore  have  two  entry  points:  an  initial  entry 
that  allocates  the  memory,  and  another  entry  that  bypasses 
these  steps.  This  is  necessary  since  the  calculator  does 
not  allow  an  array  to  be  dimensioned  twice  or  a buffer  to 
be  declared  twice.  The  initialization  routine  is  one  of 
these  two-entry  routines. 

The  selection  state  (Figs.  A-8,9)  has  a vector  to 
allow  control  to  be  transferred  to  the  TAPE  MONITOR.  The 

other  states  all  return  control  to  the  selection  state  in 

DATA  MANIPULATOR.  Each  state  will  be  considered  separately. 

TAPE  READ.  The  TAPE  READ  module  asks  for  the  name  of 
a data  file.  It  expects  the  name  of  a file  containing  data 
that  will  fit  into  an  array  A dimensioned  a[500,2].  There 
is  a possibility  of  using  this  routine  as  a subroutine,  but 

that  is  not  done  as  the  program  is  now  written  (Fig.  A-10). 

SCALE.  The  SCALE  routine  computes  a conversion  factor 
to  convert  the  output  of  the  DVM  from  volts  to  amperes. 

The  routine  interacts  with  the  operator  to  obtain  the  range 
of  the  Keithley  417  picoammeter  and  uses  the  fact  that  full 
scale  on  the  picoammeter  is  represented  by  3 volts  at  the 
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recorder  output  (Fig.  A-10) 


CONVERT.  The  CONVERT  subroutines  are  used  to  convert 
the  data  from  the  packed  format  of  the  HP-3437A  DVM  to  the 
normal  calculator  representation.  The  routines  assume  the 
DVM  is  set  to  the  10  volt  range.  If  this  is  not  true,  a 
scaling  step  must  be  used  in  the  calling  routine  to  make 
the  correction.  These  routines  are  used  in  conjunction 
with  LONG  READ,  but  can  be  used  by  other  routines  that 
operate  the  DVM  in  the  packed  mode. 

They  operate  by  dividing  the  numeric  equivalent  of  the 
ASCII  character  by  multiples  of  two,  then  using  the  frac- 
tional part  or  the  integer  part  to  form  a string  of  "r" 
variables:  rl6,  rl7,  rl8,  and  rl9.  The  sign  is  converted 

to  a + 1 which  is  assigned  to  rl5.  The  expression, 
rl5(rl6+rl7+rl8+rl9) , is  then  the  value  of  the  reading. 

LONG  READ.  LONG  READ  is  the  routine  that  programs  the 
system  DVM  and  allocates  buffers  for  the  data.  When  initi- 
ated it  samples  the  voltage  500  times.  Each  sample  is 
spaced  .05  seconds  apart.  The  routine  then  uses  the  CONVERT 
subroutines  and  a scale  step  to  get  the  readings  into  volts 
as  measured  by  the  DVM.  The  routine  then  adds  time  incre- 
ments to  the  A array.  Exit  is  to  the  SCALE  routine.  (See 
Figure  A-ll.) 

The  A array  produced  by  LONG  READ  and  SCALE  has  the  form 
of  500  ordered  pairs.  In  each  pair  (xn,yn)  xn  is  the  time 
in  seconds  from  t=0  (due  to  a problem  in  representing  In)  as 
a number,  0 is  set  as  1 x 10“6sec.),  and  yn  is  the  value  of 
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the  current  measured  at  time  xn. 

FILE.  FILE  is  the  routine  that  allows  the  operator  to 
save  the  raw  data.  It  asks  for  the  name  to  be  assigned  to 
the  file  and  checks  to  see  if  the  file  name  is  already  in 
existence.  If  the  name  is  being  used,  it  asks  if  the  old 
file  is  to  be  updated.  This  prevents  the  inadvertent  over- 
writing of  data.  When  all  these  checks  are  complete,  it 
marks  a new  file  space  to  store  the  A array  (or  uses  the 
old  file  in  the  case  of  an  update),  and  files  the  data 
(Fig.  A-12 ) . 

PLOT  LONG  READ.  PLOT  LONG  READ  produces  a plot  of  the 
data  in  the  array  A.  It  plots  both  the  I vs.  t and  the  I * t 
vs.  t curves.  For  the  I‘t  vs.  t curve  the  data  is  smoothed 
by  an  averaging  function  that  uses  the  average  of  seven  data 
points  (three  on  each  side  of  the  point  in  question)  as  the 
value  assigned  at  that  point.  PLOT  LONG  READ  does  not  change 
any  value  in  the  A array.  The  routine  scales  and  plots  the 
axis  automatically  (Fig.  A-13). 

PACK.  PACK  is  the  module  that  creats  the  B array  of 
processed  data.  This  array  is  dimensioned  B[40,3],  meaning 
it  is  an  array  of  40  ordered  triples.  The  B array  has  the 
following  form:  B(xn,yn,zn),  such  that  n has  a value  from 
1 to  40,  and  B(x^)  is  the  number  of  the  last  triple  to 
actually  contain  data,  B(y^)  is  the  area  of  the  gate  elec- 
trode  in  cm  , and  B(zi)  is  the  value  of  the  attempt  to 
escape  frequency  for  the  sample.  The  rest  of  the  array 
contains  a summary  of  the  results  of  up  to  39  uses  of 
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-12  Listing  of  FILE 
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-13  Listing  of  PLOT  LONG  READ. 


178:  "PACK": 

179:  if  Q[l,l]=0;ent  "A=  sq  cm"  ,3  [ 1 , 2 ] ; ent  "nu=  ?"  ,6  [1,  3]  ; 1+B  [ 1 , 1] 
130:  U*r6* r7 ; f or  1=201  to  500 ;A[ I,2]A[ 1,1] +r6+r6 


LONG  READ.  Each  B(xjn,yjn, zm)  where  m varies  from  2 to  B(x1) 
contains  the  following  data:  B(xm)  is  the  average  I • t value 
for  the  measurement,  B(yra)  is  the  sample  standard  deviation 
of  the  values  that  went  into  the  average  computed  for  B(xm), 
and  B (zm)  is  the  temperature  in  °K  of  the  sample  during  the 
measurement. 

After  calculating  and  assigning  the  values  as  mentioned 
above,  PACK  then  checks  to  see  if  the  B array  is  full.  If 
it  is,  PACK  asks  that  the  array  be  filed  and  exits  to  FILE  B 
If  the  B array  is  not  full  PACK  exits  to  the  selection  state 
(Fig.  A-14 ) 

FILE  B.  FILE  B does  for  the  B array  exactly  what  FILE 
does  for  the  A array,  files  it  on  tape  (Fig.  A-14). 

PRINT  B.  PRINT  B produces  a listing  of  a B array.  It 
uses  SUB  A to  find  the  array  and  put  it  into  main  memory, 
then  the  routine  outputs  the  listing  to  the  printer 
(Fig.  A-15). 

PLOT  B.  PLOT  B is  the  routine  that  plots  the  B array 
in  the  form  of  density  of  states  vs.  energy.  It  uses  SUB  A 
to  loc'te  the  array  and  put  the  array  into  main  memory,  then 
it  uses  the  equations  presented  in  the  theory  to  convert  the 
I»t  information  into  the  needed  form.  Error  bars  indicating 
one  standard  deviation  above  and  below  the  value  are  added. 
The  plot  is  automatically  scaled  and  the  axis  drawn  and 
labeled.  PLOT  B allows  more  than  one  set  of  data  to  be 
plotted  on  the  same  axis  for  easy  comparison  (Fig.  A-16). 
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rig.  a-15  Listing  of  PRINT  B. 


138:  "PLOr  B" : dsp  "PLOT  B ready" ;stp 


-16  Listing  of  PLOT  3 


T 


SURFACE  POTENTIAL  PLOT.  SURFACE  POTENTIAL  PLOT  is  the 
module  that  plots  the  change  in  surface  potential  as  a func- 
tion of  energy  by  integrating  the  charge  represented  by  the 
measured  trap  states  and  dividing  by  the  oxide  capacitance 
in  Farads  per  cm  . This  routine  uses  parts  of  the  A array 
as  temporary  storage.  The  data  is  taken  from  a B array 
located  by  SUB  A,  sorted  and  integrated  by  SUB  B,  and 
plotted  by  SUB  C.  The  scale  is  automatically  set  with  the 
help  of  SUB  S.  Axes  are  drawn  that  can  be  used  with  more 
than  one  set  of  data  points.  (Figs.  A-17,  18,  19,  20) 


Operation 

The  system  as  a whole  is  planned  to  be  self-prompting 
and  uncomplicated  to  use.  Once  the  various  elements  of  the 
experimental  apparatus  are  connected,  the  tape  cartridge 
containing  the  programs  and  data  is  inserted  into  the  cal- 


culator. The  commands,  ldp  11  i EXECUTE  1 , are  given  to 
the  system. 

The  tape  will  be  accessed,  and  the  system  program, 

TAPE  MONITOR,  will  be  loaded  and  initiated.  The  calculator 
will  respond  with: 


" Enter  the  function  code.  " 

The  list  below  provides  the  function  codes  for  the 
system  modules. 


PRINT  TABLE 

PT 

FREE  FILE 

FF 

MARK  FILE 

M 
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Fig.  A-19  Listing  of  SUB  3 


-20  Listing  of  SUB  C. 


LOCATE  FILE 


LF 


initiate  the  experiment  EXP 

continue  the  experiment  EXFC 

Type  the  code  then  press  I CONTINUE  I . 

If  the  code  selected  was  EXP,  then  the  tape  will  again 
be  accessed  and  DATA  MANIPULATOR  will  be  loaded  into  memory 
and  initiated.  The  display  will  now  be: 

" Please  select  your  option.  " 

The  list  below  provides  the  option  codes  for  the 
experiment  modules. 


vector  to  TAPE  MONITOR 

TM 

TAPE 

TAPE 

FILE 

FILE 

PRINT  A 

PRINT  A 

LONG  READ 

LR 

continue  LONG  READ 

LRC 

PLOT  LONG  READ 

PLR 

PACK 

PACK 

PLOT  B 

PLOT  B 

PRINT  B 

PRINT  B 

SURFACE  POTENTIAL  PLOT 

SPP 

In  any  of  the  routines  when  data  is  required,  the 
system  provides  a prompt.  For  questions  requiring  a yes 
or  no  type  answer  the  system  expects  a ”1"  for  a yes  and 
a "O"  for  a no. 
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The  calculator  system  and  the  software  described  here 
are  in  the  custody  of  Mr.  John  Blasingame,  AFAL/DHR, 
Wright-Patterson  Air  Force  Base,  Ohio  45433. 
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